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Summa ry . 

The following paper was submitted by the writer? as a. thecis 
to the Department of Aeronautical Engineering, at the T !arse.chusett 
Institute of Technology- It constitutes a general theoretical 
discussion of the damping factor in roll, together Tiith the re- 
sults of wind tunnel tests on the continuous rolling of a U.S.A.- 
30 airfoil- Tro general formulas are derived for +he damping of 
roll, each of which contains unavoidable indeterminate functions. 
Certain of these functions have been evaluated from the test data- 
Of chief interest is the deduction that the actual damping ar ex- 
perienced in flight differs from the damping as theoretically 
calculated ny a function of the ^ing-tip pressure distribution, 
which is in turn largely influenced by the fonri of the win^-tip 

and by the rolling velocity. Finally, it has been shora that in 

dCr d0 7 
the damping equations may he substituted for even 

under full flight conditions, .without serious error. 
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FART I- 
TREORHTICAL DISCUSSION. 

Ijltrgductf. on . 

The damping of an airplane in roll is an exceedingly coumlex 
problem and one which it may never be possible to solve completely, 
but if 7/e can form some idea of the magnitude of the principal 
factors involved and of their importance in practical flying, a 
mathematically complete solution can be dispensed with. Our int- 
erest in the damping coefficient is mainly as a guide to the de- 
termination of stability and as an indication of the forces en- 
countered in maneuvers. 

The damping coefficient itself is made up of many components, 
arising from the several elements concerned: v:ings, fuselage, tail, 
etc. As long, however, as we confine our tests to the complete 
airplane, T .ve have no means of analyzing the source of the magni- 
tude of the basic elements of damping. The chief contributor to 
these elements, the r ^ing, is itself affected by a series of com- 
plications. Dihedral, stagger, taper ?nd aspect ratio are only a 
fev; of the complicating factors. In view of the foregoing, there- 
fore, this report -ill be confined principally to the straight 
rectangular wing of constant chord and constant section. 

Theory of Dimens ions . 

Probably the simplest method of attacking the problem is by 

the theory of dimensions, ^e may express the damping coefficient 

of a rectangular wing of constant section as a function of several 
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variables, thus: 

,dCz . v . . 

iJT = p = f ficT> u > p ' c > ^ P> 

I> = the total damping coefficient 
arising from all causes talc en 
together, as distinguished from 
L p , the damping coefficient of 
roll due to roll. 

L = the total rolling torque in lb-ft« 

p = angular velocity of roll in radians 
per second. 

dC 7 

= rate of change of the coefficient of 
normal force with angle of attack 
(i.e., the slope of the normal force 
curve) . 

The units of da are radians. 

u = air speed in feet per second. 

c = wing chord in feet, 

b = wing span in feet. 

z = perpendicular distance from the 
axis of rotation to the mid- 
point of the wing chord, in feet* 

P - mass density of the air. 
dC 7 

The term -^a is used instead of the angle of attack because 
the same angle of attack does not give the sa;;<e results for differ- 
ent wing sections, or for the same wing section tested at differ- 
ent values of VL. Thus we eliminate both scale effect and the 

dC 

effect of different wing contours- is the slope of the curve 

of that component of force perpendicular to the line of steady 

flight, (the trajectory of the center of gravity) and should not 

dCr 

be confused with ■ which is the slope of the lift curve as 
read from a wind, tunnel plot and relates to the force perpendicu- 
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lar to the relative wind, whatever direction that may take* 

It is ucually assumed that the torque caused by a force on 
an elemental area do, is proportional to C z d3 or to 
(dCz/da) (Aa)dS, and it would therefore seem necessary to intro- 
duce a Ac term into the foregoing equation. Aa is the differ- 
ence between the real angle of attack at the point in question 
and that at the plane of symmetry. However, since Ac is always 
small, Aa^ tani£a) = py/u, y being the distance from the plane 
of symmetry along the span, and since both p and u are al- 
ready expressed in the equation it not be necessary to add 
any A a term* 

Evaluating the various elements in the above equation by the 
theory of dimensions, v;e get: 



Lrp = k 



dC * , 4 
a u 



f ft» (°\ f-M (g\ (1. 

V U / \C/ yby 



The first of these indeterminate functions contains the 
ratio of the linear velocity of rotation of the ring-tip to the 
wind velocity which amounts to a particular value of tan(^a) 
or act, approximately. The second contains the aspect ratio. 
The third contains the ratio of the height of the wing above the 
rotational axis to the span of the wing. Note that all of these 
functions are dimensionless ratios. 

Before we can make any use of this equation it will be nec- 
essary to examine the three indeterminate functions more fully, 
and check the results of our theory with the experimental data. 
To this end a more elaborate, if less convenient theory has been 
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developed • 
General Theory. 

The damping coefficient in roll is usually calculated theo- 
retically by assuming the wing to be made up of weny minute ele- 
ments, the lift on each of which contributes an element of rolling 
moment. Thus in the accompanying sketch, 
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L* = the theoretical rolling torque- 

dS = element of area = chord x dy = c(dy) ■ 

V = resultant wind velocity where 

V 2 = u 2 + (pR)' or V 2 - u 2 + p 2 y 2 + p 2 z 2 

cl c 

C 2 = normal force coefficient = ( M 

where g^ 21 is the slope of the normal 
force curve, and is assumed constant, 
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while Ac in the chanpe of ar^le of 
attack along the win? from the plane 
of ^vnTnetry. is the cliff erenco 

between the value of C z at the viaz 
center line and that at any point, 
y along the span* 



Therefore: , • 



T ! — O 

- - • d 
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At any point y elong the wing A a tan - py/u, and since 
Aft is necessarily small we may take £ci - py /u, where Ac is 
in radians, without appreciable error. 



Then: 



Mb/ 2 

/ (£2L ay \ (u £ + p 8 y 2 + p 2 z*)dy 

Jo da ^ u ' 

Integrating, and collecting terms, we obtain 



L- - (pcb 3 ; 



(„ + Sill) + ;_L (B^-) 



ia u ' 80 u 

This formula neglects the irregular pressure distribution 
at the wing-tip, or any change in that distribution due to the 
rolling motion. If we generalize the equation of roll so n s to 
include the tip effects and to include the effects of the other 
rotary and resistance derivatives as well, we get the total roll- 
ing torque L, thus 

L = L« + [f(t)] b + Lv v + Y v v z + Y p p z (3. 

where L 1 is the theoretical rolling moment due to roll obtained 
by equation #2; f(t) is an indeterminate function of the tip pres- 
sure distribution, and produces roll by acting in the direction of 
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lift, at a distance from the axis of rotation vrbieh is some frrc- 
tion of t; : e S'oar b; L v is the rolling ncmenc duj to side-slip, 
therein the side-slip is introduced by the fact that the ?:ing ic 
rotating at a normal distance s fro in the axis of roll: v is 
the amount of the side-slip velocity; and Y~ and Y*> are the 
lateral forces due to side-slip and to roll respectively, which 
produce roll by acting along the ning span at an arm z with 
respect to the axis of rotation. 

If r;e substitute v = r>z in equation tt? (where p is ir 
radians per unit time) and divide through by p, wc obtein 

p p J v " * 

and since torque divided by the corresponding velocity gives the 
coefficient, we have the general coefficient of roll depressed 
thus: 

L T = L p + [f(t) p ] b + [L v + Y p ] ? + Y v z* (4. 

In this equation the fCt^, serves as a correction factor for 
the theoretical damping of roll (L'p' , which night be expected 
to be too large since it neglects the falling off of the lift at 
wing-tips- The terms L v , Y v> and Y^ are the experimental val- 
ues for the wing in question- If no-? we divide equation #2 by 
p, we get an expression for L'/p fi.-e., L*p) which may be sub- 
stituted into equation #3. Collecting teri^s, v/e have the com- 
pletely general equation: 
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L T = Tl (da , cb u + • 80"Ua ) u 

[f(t) p J b ♦ [L v + Y p ] z + [Y v + ^ z»] (5. 

clC 2 

In this equation the three terms involving jjo"- and the 
term [f(t) p ] b together make up what is usually called L p , 
damping of roll due to roll. The first of these 1 % ^^cb 3 u, 
represents the damping obtained by the element theory, if the re- 
sultant airspeed is everywhere taken equal to the speed of flight- 
The second term represents the added damping obtained if wc tn\e 
the airspeed as the resultant of the speed of flight and the nor- 
mal sneed due to rotation, and, under full flight conditions 
with the maximum probable velocity of roll, is about 2$ of the 
first term. The last term -^(|^)p 2 cb 3 z 2 represents the 
further increment of roll added by considering the transverse 
component of airspeed across the wing, and amounts at most to ap- 
proximately 0.25$ of the first term- The other terms have al- 
ready been discussed. 

Discussion of the General Theory . 

Obviously the general equation is too complicated for conven- 
ience and will have to be simplified by assumption. If we neg- 

dC 

lect the two smaller terms in jo^ and tnereb y introduce an error 
of not more than 2 l/2$, the equation becomes: 

L T = T2 (1lf) cb3u + Cf(t) p 3 b + 0* + Yp] z * Y v z 2 (6. 
Further simplifications will depend upon the specific conditions. 



If 2 is small, "both of the last terms disappear; if z is ma- 
terial, Y v is usually negligible; Lv and Yp increase with 
dihedral, hut the latter is apt to be unimportant. In the case 
of the rectangular wing of constant section, all three terms 
L v , Y v and Yp, can probably be neglected and the equation takes 
the form 

l t B Ire* cb * u + tf(t) p ] b (?. 

It will be noted in this case that Lj = L p since the correction 

factor [f(t) p ] b = L p - L' p . 

A study of equations Nos- 5, 6 and 7, will shed some light 

on the variation of L with changes in aspect ratio and wing area 

dC 

In equation #5 the first term and the last term involving de- 
pend upon cb 3 , while the second term depends on cb 5 . The third 
term [f(t) p j b, is rather difficult to analyze. Being a func- 
tion of the tip pressure distribution, f(t) p evidently depends 
upon the chord; also the extent to which this irregular distribu- 
tion extends inward will presumably be a function of the chord 
rather than of the span- We might say, then, that [f(t) p ] b 
depends on c 2 b. However, the question arises as to whether the 
form of the tip distribution does not depend upon the normal com- 
ponent of the wing-tip velocity, which, in turn depends directly 
on the span. In other words, does not [f(t) p ] b depend primar- 
ily upon c 3 *b s ? The latter seems more reasonable. Of the remain- 
ing terms, Lv depends upon the area, the amount of the dihedral, 
and the span, or on cb 2 . Yp and Y v do not depend upon the 



span to any appreciable extent, but rather on c 2 . Summing up, 

therefore, it appears that the most important term depends on cb 

the term [f(t) p ] b, which may amount to r 60% of the total depends 

1 /uC... N 

either on c s b or c 2 b 2 ; the term involving g -Q yfix / » vmich 
comprises about 2$ of the total, depends on cb 5 ; and the almost 
negligible terms Yp and Y v depend on c 2 . 
Applying the formula 

1 



b^ = 



(be) J 



x+y 

i 2 



y_-x 



where S and R are area and aspect ratio respectively, it fol- 
lows that the cb 3 terms depend on S 2 R; the c 2 ^ term on S 2 
alone; the c 2 b, if re choose to use it, on G 3/2 R~ 1/s ; the cb 5 
term on S 3 R 2 ; and the two c 2 terms on SR" . At first eight 
the term involving S 3 R 2 would assume undue importance. Actu- 
ally this term also involves p 2 (the square of the rolling ve- 
locity) which obviously decreases at about the same rate, ox, per- 
haps, faster, than S 3 increases in actual flight, so the propor- 
tion of Sfo of L<p which was obtained for that term for typical 
flight conditions on a 2000-pound airplane, will probably not be 
exceeded for airplanes of any size or proportion. 

Neglecting these less important terms, we come to a study of 

2 

equation #7, wherein the first term depends upon S R and the 
second upon S 2 or S 3/S R -1/2 depending on how we consider f(t) p 
It seems most reasonable to take [f(t) p ] b as dependent upon 
c a b a , (or on S 2 ) which has the added advantage of bringing in 

4 

the area without fractional exponents, and checks the b term 
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obtained "by the theory of dimensions in equation #1. It has been 
found by the N.A.C.A.,* however, that in going from model test to 
full flight, Lt increases less rapidly than S 8 , so there is 
something to be said for the alternative supposition. In either 
event, the mean exponent of the aspect ratio is bound to be less 
than unity, probably around .8, since about 70$ of L T depends 
upon aspect ratio to the first power. This also has a bearing on 
equation #1. 

The Slope of the Normal Forc e Curve- 

dC z 

Throughout the discussion we have used the term ratner 
than g^k. If we take lift as perpendicular to the relative wind 
and the "z" force as perpendicular to the line of flight, we have 
by the familiar transition 

C z = Oleosa + C D sino 

Differentiating, 

= 5^cos(Aa)-CLsin(Aa)+C D cos(.Aa)+3-d n sin(Aa) 

• (8. 

(aa) here represents the change in angle of attack from the value 

at the wing center-line- ( Act) is zero at the center-line. At 

dC z dC L 

that point, therefore, gcr = cTcr + C D- 

A full-scale example has been worked out in Fig. 1, for a 
U.S. A. -30 wing of 60 ft- span and 10 ft. chord, turning at 1-5 
radians per second, which is certainly an exaggerated case- It 
will be seen that the 3^ curve follows the curve very 

* N.A.C.A. Report #167 
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closely. Above a = -6° the deviation is less than 3 \/?jfo througn- 

out the entire range. This means that for all present-day air- 

dCr, dO Tl 

planes under ordinary conditions d'a"" - ST'* ?or auto-rotation 
the more exact form is required, and a new curve of normal force 
must be plotted for each angle of attack considered, and a graphi- 
cal solution must be made if the roil is very rapid. 

Experimental Results * 

As a conclusion to the theoretical discussion, the results of 
the experiments may be summarized- The indeterminate functions 
of {^) and (|u0 in e Q uai:ion w e;re investigated for an aspect 
ratio of 6. It was found, as might be expected, for the plain 
rectangular wing of constant section, that L-j was independent of 
(^yt at least within the experimental error- This follows from 
equation #6, since L v , Y p , and Y v are known to be small. The 
other indeterminate ^(^f*^ takes the form of k T (^p) > where 

k a and n vary practically as straight line functions with 

dC 7 / dC 7\ * /r>b\ n 

j^^. The equation then takes the form Lt = kj (cfo j ub * 

Values for k a and n are plotted in Fig. 2. The full lines rep- 
resent decreasing values of and cover the range of an aver- 

age lift curve from maximum steepness up to nearly maximum lift. 

dC z 

The dotted lines represent decreasing values of which cover 

the range of the lift curve below maximum steepness (i-e., in gen- 
eral, below a = 0°) where the lift curve tends to bend upwards 

from a straight line. The point of intersection is at the maximum 
dC 

value of jj^ 2 -* Obviously, with wings having a different maximum 
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vaxue of > *^ ne intersection would take place at some other 
point, and presumably the dotted line would be displaced vertical- 
ly to correspond. The fact that these lines intersect at n - 0, 
for this particular wing, is probably accidental. However, the 
full lines, which cover the normal range of flight angles, should 
apply equally well to any plain rectangular wing of constant sec- 
tion and aspect ratio 6. The curves show very clearly a depart- 
ure from a linear relation between rolling moment and rolling 
velocity when.: the rolling velocity is high. 

A value of can be calculated from the element theory on 

the assumption that n is equal to 0. Using equation #2 and neg- 
lecting all terms within the bracket except the first, we have: 

1 /dQ„\ 
L P - T2~ (da / c ^ 

where the units are homogeneous throughout. In plotting the curves 

dC 

of Fig. 2, ^— was taken in Ib/sq. ft/mile per hour/degree, and 

u in M.P.H. This introduces a correction factor of 57.3 x 15/22. 

A further factor of 1/6 is introduced by the substitution of b* 

for cb 3 . The total value of the calculated coefficient, on a 

basis comparable with that used in Fig. 2, is therefore, 57.3 x 

15/22 x 1/6 x 1/13, or .52. The experimental constant will be 

seen to approach the calculated one closely at small values of 
dC 

^■g 2 -, but it falls far below at high values of the slope (corre- 
sponding to angles of attack well below that of maximum lift). 
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PART II. 

EXPERIMENTAL DETERMINATION OE THE DAMPING OF ROLL. 
Metho d of Tes t . 

The tests were carried out in the 4-foot wind tunnel at the 
Massachusetts Institute of Technology/, on an 18 :! x 3 tJ wood model 
of the U.S. A. -30 wing at a wind speed of 40 M.P.H. A 13" spindle 
1/2" in diameter was mounted axially of the tunnel between a pair 
of conical bearings. Each bearing was supported by three wires 
to the side of the tunnel, so attached as to keep the bearings 
seated snugly on the ends of the spindle. The spindle was provid- 
ed with a pair of slots 2 V apart, through which passed two l/8 u 
diameter rod3, the rods being screwed into the wing model at mid- 
span, one behind the other- Special counter-weights were mounted 
at the opposite ends of the rods. Thus, by loosening two set- 
screws opposite the spindle slots, both the angle of attack of the 
wing and its distance from the axis of rotation could be altered. 
Finally, a light flexible cord was wrapped three times around the 
spindle and the ends carried out through the bottom of the tunnel 
to a pair of weights, which supplied the driving torque to keep 
the model in continuous rotation- The mounting is shown in Fig. 6. 

Runs were made at various torques for five angles of attack 
(-4.5°, 0°, 6°, 12°, and 18°) and at 0° for four positions of the 
wing relative to the axis. Each run was reocated in the reverse 
direction and the results averaged to remove any error due to warp 
in the wing. The speed of rotation was observed for each \ r alue of 
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torque by counting the revolutions of the wing against a stop-watch. 

The net toroue was then obtained by subtracting a friction and 

windage correction which has been found independently by experiment. 

Torque in foot-lb divided by angular velocity in radians per second 

gave the damping coefficient. Tests were also made to determine 

the usual characteristics of the wing and are given in Fig. 3. 

dC 

From these tests the values of t— 2 - were obtained. 

da 

Analysis of Results . 

In order to determine the unknown function of z/b in the di- 
mensional equation #1, values of L T were plotted against z/b in 
Fig. 4 for three values of pb/u. Unfortunately, it was not possi- 
ble to make pb/u exactly constant, since in making the experi- 
ment the speed of rotation for a given torque could not be fore- 
told. However, groups of values were selected in which pb/u is 
essentially constant and since all values for all groups lay with- 
in 5fo of the average, which is within the error of the experiment, 
the evidence seemed sufficient to indicate that L is practically 
independent of z/b for straight rectangular airfoils. A few 

tests were made at a = 12° for various values of z/b to deter- 

dC 

mine the effect of a change in ^g 2 - on z/b. These results were 
slightly more erratic than those at 0°, but nevertheless bore out 
the fact that is independent of z/b- This is what the theory 

led us to expect, knowing that L v » Y p , and Y v are small. 

To evaluate the pb/u term in equation #1, the runs for each 
angle of attack were plotted on logarithmic paper, using pb/u as 
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abscissae and LT/(dCz/da)ub* as ordinates- It was found that 
e r ..cl'i :-e ij of . mints lay nearly along a stiaight line. Accordingly, 
from the intercepts of this line the values of k a and n "riven in 
Fig- -3, were obtained. These values \ieie plotted egainrrt the 
slope of the lift curve instead of the angle of attach so as to 
be of more general application* 

It is interesting to compare the theoretical and experimental 
values of L. Fig. 5 shors a typical comparison. The hotter, line 
represents the experimental values of L<j and the top line the 
theoretical values of the conventional Lp obtained from equation 
#3- Since both theory end experiment agree that for the straight 
rectangular wing the extraneous terms Lv, , etc, are negligi- 
ble, it follows that the only difference between the theoretical 
L p and the actual L T should be the tern [f(t) p ] b, (i.e., the 
tip loss correction) . With this in mind the central line in Fig. 
5 was obtained by solving graphically for Lp and assuming the 
conventional tip pressure distribution of l/3 the running load at 
the tip, tapering up to full load at .5 of a chord-length inboard 
from the tip. Evidently, then, the tip pressure distribution is 
altered by the rolling motion. This is not unexpected. 

The comparison shown in Fig. 5 gives the greatest deviation 
Tvhich was found. As the angle of attack is increased, the devi- 
ation becomes less, until near maximum lift the experimental value 
becomes the greater of the two, as already noted in Part l } in 
connection with the discussion of Fig- 2. Except for the varying 
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distortion of the tip loading no adequate explanation accounts 
for this peculiarity. 

Conclusi on. 

Finally, it must be recalled that the foregoing is in the 
nature of a theoretical discussion, and that the experimental data 
represents only a single wins: on which we cannot afford to gener- 
alize too much. We must have further data. Specifically, we re- 
quire wind tunnel tests on wings with dihedral and taper, on bi- 
plane combinations, and on different tip forms; we require free 
flight tests on various airplanes, both large and small, especi- 
ally monoplanes. T/ithout these additional data very little tangi- 
ble progress can be made. 
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Full lines for decreasing values of dC-/do. Dotted lines 
for increasing values. dC_ in lb/sq.f t/M.P.H. Units 
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Theoretical values of rolling due to roll 
(L, ), obtained from equation No. ° 




Theoretical values of rolling due to roll 
found graphically with allowance made for 
tip loss. 



Experimental values of total rolling factor 
(Lj) obtained from wind tunnel, j 



,06 



.OS 



.10 .12 

Values of pb/u 



.16 



dC z /dc£.000lSfc lb/ sq.ft/M.P.H./degree=. 00^91 lb/sq.ft/ 
f >p. s. /radian 

Fig. 5 Comparison of theoretical and experimental values of 
the damping factor in roll. 
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